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Abstract. Surface codes describe quantum memory stored as a global property 
of interacting spins on a surface. The state space is fixed by a complete set 
of quasi-local stabilizer operators and the code dimension depends on the first 
homology group of the surface complex. These code states can be actively 
stabilized by measurements or, alternatively, can be prepared by cooling to the 
ground subspace of a quasi-local spin Hamiltonian. In the case of spin-1/2 (qubit) 
lattices, such ground states have been proposed as topologically protected memory 
for qubits. We extend these constructions to lattices or more generally cell 
complexes with qudits, either of prime level or of level S for d prime and f > 0, and 
therefore under tensor decomposition, to arbitrary finite levels. The Hamiltonian 
describes an exact Z^j = Z/cZZ gauge theory whose excitations correspond to 
abelian anyons. We provide protocols for qudit storage and retrieval and propose 
an interferometric verification of topological order by measuring quasi-particle 
statistics. 
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1. Introduction 

There is a rich history to the study of topologicaUy ordered states of matter. Such 
states are defined by the property that aU physical correlation functions are topological 
invariants. In the field of condensed matter, these states have been proposed as ground 
states of models for high temperature superconductors and for fractional quantum Hall 
states [22] . Furthermore, it has been demonstrated that such order can arise as a low 
energy property of hard core bosonic spin lattice models. In contrast to the familiar 
situation with spontaneous symmetry breaking, here the ground states exhibit more 
symmetry than the microscopic equations of motion. It has been suggested that such 
emergent properties may model gauge fields and particles found in nature [16]. In 
the field of quantum information it was shown by Kitaev [12] that ground states of 
Hamiltonians which can be expressed as a sum over quasi-local stabilizer operators 
provide for topologicaUy protected qubit memories. These states are referred to 
as surface codes. They are robust to arbitrary quasi-local perturbations and have 
abelian anyonic excitations. In order to perform universal fault- tolerant quantum 
processing, it is necessary to use non abelian anyonic excitations that transform under 
an appropriate group [12]. From the algorithmic point of view attempts have also been 
made to understand quantum computing in terms of nonabelian anyon operations 
[1]. Finding suitable microscopic lattice models that provide for universal quantum 
computation is an area of active research [(!]. 

This survey attempts to exhaust the topic of surface codes for topologicaUy 
protected qudit memories. While not as powerful as fault tolerant models with non 
abelian anyons, these models offer a new perspective on non-local encoding of quantum 
information and give us insight into microscopic realizations of lattice gauge theories. 
Surface codes for two level systems [12] are by now well understood. Their implications 
for error-resistant quantum computer memories have also been considered [,]. In the 
error-correction context, the topologicaUy ordered eigenstates may be understood as 
a particular case of quantum stabilizer codes (e.g. [9].) The error lengths of the 
resulting stabilizer codes are not exceptional, and only rarely do anyonic systems 
appear in classifications of near-optimal quantum codes. (Optimality in this sense 
refers to minimizing the number of code-qubits against the number of errors a code 
may correct.) Yet all the error correction operations are local upon the lattice in 
which the quantum data is stored, which might improve scalability. Moreover, an 
aside to an argument focused on deriving a famous stabilizer code from the topology 
of the real projective space in fact demonstrates that a qubit lattice is not required [s]. 
Rather, a two-complex (see e.g. [15]) suffices, where a two-complex is a generalization 
of a graph in which discs are also allowed with edge boundaries. On the physical 
system which places a qubit on each edge of a (cellular or simplicial) two-complex F, 
there exists a Hamiltonian whose topologicaUy ordered (stabilizer-code) groundstates 
are parametrized by the first homology group of the complex with bit-coefficients: 
Hi{T, F2). The Hamiltonian is a sum of vertex and edge terms which are proportional 
to cither tensors of Pauli Z operators around qubits on edges adjacent to a vertex or 
are proportional to tensors of X operators on edges bounding a face of the complex. 

For some time the existence of stabilizer codes over qudits (d prime) have been 
known [•)]. Yet only recently have results on the topic become as strong as those 
applicable in the bit case, including estimates of optimal code-lengths etc. ([11], see 
also [10].) Moreover, extensions to prime-power (d^) level qudits (actually qud^its) 
have also been found, so that tensors provide a stabilizer formalism for all finite-level 
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systems. In this work, wc exploit the new stabiUzer formaUsm to construct codes on a 
two-complex whose edges carry prime-d-level qudits, and we also outline the extension 
to d^-level qudits. The associated ground states are parametrized by Hi{r,¥j^e), where 
the coefficient field is viewed as an abelian group under addition. This requires few 
new ideas, although care must be taken with sign conventions which were vacuous in 
the earlier work on F2-coefficients. Thus, after tensoring we have constructed surface 
codes with qudits for arbitrary finite d placed on the edges of a generic orientable 
two-complex F. Recent work by Bombin and Martin-Delgado [2] investigates classical 
and quantum homological error correction codes. They construct a class of surface 
codes for qudits which asymptotically saturates the maximum coding rate and provide 
several example encodings on various two complexes. Here we do not address the issue 
of coding efficiency. Rather we concentrate on explicit constructions of Hamiltonians 
that support qudit surface codes in their ground eigenstates and describe how one 
might encode and decode therein. 

The manuscript is intended to be self-contained. Thus, §2 opens by reviewing 
some the required facts on stabilizer codes. In order to aid readers less interested in 
the general case, §3 treats prime-d level encoding on surfaces separately. Methods for 
encoding; decoding, and stabilizer measurements are given in §4. Extensions to the 
case of prime power qudit encodings arc given in §5. Errors in our model correspond to 
low lying excitations in the Hamiltonian whose superselection sectors may be viewed 
as massive particles on the underlying cellulation. In §6 it is shown that our model 
reproduces a Z/dZ gauge theory where errors are described by particle anti-particle 
pairs of change/flux dyons. We propose an interferometer circuit for measuring the 
the statistics of these quasiparticles. We conclude with a summary and some. 

2. Qudit Stabilizer Codes 

We next review stabilizer codes [9, 11]. This section focuses on the case of qudits 
with a prime number of levels. The first subsection recalls the definition and a basic 
technique. The next subsection generalizes a well known construction from bits to 
dits. 

2.1. Stabilizers and groundstates 

Let d be a prime number, and consider the qudit state space 7i(l, d) = C |0) • • • 
C|d— 1), with a pure state of n qubits being a ket within H{n,d) = 7i(l,d)®". A 
possible generalization of the Pauli operators on 7i(l,d) would be to consider the 
group generated by the following unitary matrices: 



These are not Hermitian unless d = 2. The qudit Pauli-tensor group, say V{n,d) C 
U[H(n,d)], is the group of unitary matrices generated by ri-fold tensors of elements 
oi{Id,X, Z}. 

Wc might be more explicit in the description of 'P{n, d). First, for n = 1, label the 
multiplication in to be a dot-product. Then Z'^X" = ^'^•''X'^Z'' . More generally, 
for dit-strings a, be (F^)", we use X®"^ and Z®'' to abbreviate X"^ (g) X"^ (g) • • ■(g)X''" 
and similarly Z^^ for Z^'^ (X) Z^^ ■ ■ ■ Z^" . For the n-entry dot-product with values 
in Fd, we have Z'^^'X'^'' = ^<^*b x®°- Z^^ . Thus explicitly 



X\]) 
Z\3) 




for ^ = exp(27ri/d) 



(1) 



V{n,d) 



Z®"; a,6e (Fd)",ceFd} 



(2) 
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The qudit stabilizer groups are subgroups G C 'P{n,d). The code subspace of such a 
stabihzer group is the joint +1 eigenspace of all g & G. 

Of course, such joint eigenspaces might well be trivial. Yet a standard argument 
shows that they are nontrivial in certain cases. This technique is so fundamental 
to stabilizer code manipulation that we wish to highlight it; it will be used several 
more times in the course of the work. While actually an elementary technique from 
representation theory, it has also featured prominently in the quantum computing 
literature [1 4]. 

Stabilizer code projectors: The sum of unitary maps tt = (#G)^^ Sgec^ ''^ 
projection onto the code- sub space. We present the argument. First, tt^ is the identity 
map since irg = tt for any g E G. Second, tt = tt^ since adjoints are inverses in 
the unitary group. Hence, either tt is a projection or —1 is an eigenvalue of tt. 
Yet Id" is a summand, so the complex inner product precludes ttIiP) = — {ip) for 
a nonzero Now split T-l{n,d) = Vi © V2 © • • • © 14 into irreducible orthogonal 
unitary subrepresentations of G. For each Vj, the image under tt and its orthogonal 
complement form a decomposition of Vj . Thus by irreducibility, tt either preserves a 
Vj or TTVj = 0. Clearly the former holds for any irrep (i.e. irreducible representation) 
within the code subspace of G. On the other hand, if {'ip\g\4') 7^ 1 for some g, then 
the latter holds. 

As a remark, irreps within the code subspace of G must be one-dimensional and 
are also known as trivial representations. As a second remark, the code subspace is 
nonzero iff Trace(7r) 7^ iff (G n {^Ud}) = {h^}. 

In the Hermitian case {d = 2,) it is standard that all eigenvalues of group elements 
are ±1, so that a suitable Hamiltonian for which the code space is the groundstate 
is — TT. For general d, the eigenvalues lie within the unit circle, so that —1 is still the 
least possible real part. Also, |A) = (1/A) |A) = A|A) since ~ g~^ . Thus, one 
may place the qudit code subspace into the groundstate of a Hamiltonian by adjusting 
each summand of tt with a Hermitian conjugate: H = X^gec "I" 5^)' that the 
eigenvalues of the summands are then — 2Re[spec((7)]. 



2.2. Quantum circuits for qudit stabilizer measurements 

Given an n-qudit system, it is important for purposes of error-correction to be 
able to test whether or not a state j-f/;) lies within the stabilizer code of some 
^ = ({.9i}) ^ 'P{n,d). It suffices to test whether is a -1-1 eigenvector of each 
generator gj. We sketch quantum circuits which achieve such a measurement. 

Let J-d ~ d'~^/'^Yl'j~h=oi^'^ b) (^1 be the qudit Fourier transform. Considering 
eigenkets, T\XTd ~ Z. Now the number operator n = X]^=oi li) 01 suffices to infer 
the eigenvalue of Z and project into the appropriate eigenstate. As a circuit, we might 
denote a number operator measurement with the Z symbol, one of several common 
conventions in the qubit case: 



Determination of the X eigenstate may be accomplished by 



Similarly, there is some one-qudit unitary which will diagonaHze any X°'Z^ € 7^(1, d), 
usually not a Fourier transform. Yet using the diagonalization and a number operator 
one may infer an eigenstate. 
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For Z®*^ and X®*^, wc suggest using addition gates along with a qudit ancilla. 
We will denote \j, k) ^ |j, (j + k) mod d) by a typical control bullet with the target 
(in the formula second) line holding a + gate. The the following construction of Z®^ 
generalizes for Z®^: 

|0) ■ 



For Z'^^ \j1J2, . . . , jfc) = ^-Ji+ ' +J" \j1J2, ■ ■ • , jfe), and we have placed |(ji H h jfe) mod d) 

on the ancilla line before the number operator is applied. Note that Z ® Z~^ results 
by replacing one of the modular addition gates above with modular subtraction. Pow- 
ers of operators are measured by multiple appications of the sum gate appropriately. 
Finally, (jrt)®fcj^»fe_;r»fe ^ 2■®^ so that the following diagram for extends: 

|o) — [+]-{+]-{^= 



Td - 



Using similarity transforms by qudit Fourier transforms, we may similarly achieve 
X ® Z ® Z ® X etc. Yet more generally, the comment on existence of diagonalizations 
above produces circuits for arbitrary elements g V{n,d). 



3. Homologically Ordered Groundstates for Prime Qudits 

It is typical to place topological orders on explicit planar or spacial lattices of spin- 
j particles, e.g. square, triangular, hexagonal, Kagome, etc. An alternative was 
presented in Freedman and Meyer's derivation of certain error-correcting codes of 
Shor and LaFlamme [s]. Namely, qubits could be placed on the edges of a two- 
complex F, and an appropriate Hamiltonian would have the dimension of its degenerate 
groundstate eigenspace equal to the number of classes within Hi{T,¥2)- We next 
extend this construction to prime- level qudits; the task is mainly to keep track of sign 
conventions which are vacuous in F2. We then check that the groundstate eigenspace 
is similarly spanned by kets associated to elements of i7i(r, F^), by applying stabilizer- 
code techniques. 



3.1. Cellular Hamiltonians 

Label V to be the vertices of F, £ to be the edges, and !F to be the faces. We also 
require properties that hold if F is a cellulation of an orientable, compact, connected 
surface. Specifically, each edge has a boundary of exactly two vertices and each face 
has an orientation according to which each edge lies in the boundary of two faces 
with the edge taking opposite orientations in the boundary of each face. Finally, F 
is finite and H2(T, F^) is a copy of F^ spanned by [F], the sum of all faces with their 
orientation according to F. 
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We briefly review the appropriate homology. Label the chain sets to be formal 
sums of vertices, edges, and faces respectively: Co(r,Frf) = span][r^(V), Ci(r,¥d,) = 
span]f.^(f), and Ci^^^d) = span^^ ( J^) . We generally drop the F and coefficient 
system, which should be clear from context. Since F is a cell complex, there exist 
boundary operators 

Co^ Ci ^ C2 (3) 

with = [I'l]. For example, if an edge e connects vi and V2, say e = [wi,f2], then 
de = vi — V2 = vi + {d — l)v2- Note that the definition of F demands that edges 
e G £ are images of [0, 1] within F, and hence all edges are implicitly oriented. The 
coefficients further allow for F^-valued multiplicities on each edge. Since = 0, we 
have ker(9i) 3 image(92) for dj : Cj — > Cj-i. Thus we may define the F^j vector 
space Hi{T,¥d,) ^ kcr(9i)/image(92). This first homology group is well known to be 
a topological invariant, i.e. any topological space homotopic to that underlying F will 
produce an i?i(F,Fd) of the same dimension. Homology elements are represented by 
cycles, i.e. elements of the kernel of the boundary operator. However, several elements 
might represent the same class, differing by a boundary, i.e. an element of d{C2)- 

Recall that any Hamiltonian on n-qudits may be written as a sum of tensor 
products of Hamiltonians (Hermitian matrices) on each factor. The degree of a 
summand in the tensor basis is the greatest number of non-identity factors in any 
term. A fc-local Hamiltonian is a Hamiltonian whose degree is bounded by k in some 
decomposition. The topologically ordered Hamiltonians defined below are fc-local for 
k the maximum of the valence of any vertex and the number of edges on any face. 

Let n = and consider placing a qudit on each e G £. Again, each edge is 
the image of [0, 1] and is oriented (by F) from one vertex to the other. For the qudits 
associated with each edge, the |1) excitation of the edge will be implicitly associated 
to this orientation, while the \d— 1) state corresponds to the other. 

On the associated physical system Ti,{n,d), let and denote the operator 
applied to the qudit of that edge with identity operators buffered into the remainder 
of the tensor. For each u e V, we define a Pauli-tensor and vertex Hamiltonian by 

9v = Y[e=[*,v] ne=['U.*] ^ / 4\ 

= ~{9v + gl) ^ ^ 

For some [/ > 0, we then define the potential energy term of a topologically ordered 
Hamiltonian by Hq = VYliv^v 

The notation Hq has been chosen for the following reason. Suppose that uj = 
Seef ""eG is a chain, with each G F^. There is an associated qudit computational 
basis state, say jw), which is local and places the qudit of each e in state jrie). We 
claim that |w) is a groundstate of Hq iff dui = 0, i.e. w is a cycle. To see this, one 
verifies that gu |w) = (.'^ \uj) where duj = cv + J^w^v CwW. Hence jw) is in the stabilizer 
({ffti}) ^ 'Pi'n, d) iff \uj) is an eigenstate of each Hy of minimial (real) eigenvalue iff \oj) 
is in the degenerate groundstate eigenspace of Hg. 

Strictly speaking, one should not refer to the groundstate of Hg as being 
topologically ordered. Admittedly, groundstates are of the form \tpg) = ^ a„ 
for UJ a cycle, colloquially a loop of excited edges. For d > 2, the edges must be 
properly oriented, and hitting every edge of a y junction is allowed if multiplicities 
are accounted for. Yet the cycle subspace is not a topological invariant. Indeed, should 
F be a cell complex, subdividing F by breaking each 2-simplex (triangle) into several 
subtriangles will generally increase the size of ker(9i), although such a subdivision 
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Figure 1. Cellulation of an oriontable surface. Each system particle (qudit) is 
represented by an edge. Particle interactions occur between all edges that meet 
at a common vertex and all edges comprising a plaquottc boundary, (a) In this 
example, physical qudits reside on the vertices of a Kagome' lattice on a torus such 
that the resultant cellulation is a honeycomb lattice on a torus. Edge and face 
orientations are indicated. For the vertices vq,v\ and faces fo,fl the mutually 
commuting operators in the Hamiltonian are g„p = Z[„g ^^j, = 



,X7 



^X\, 



9fi 



Xt, 



(b) Same cellulation with 



vertex (red) ancilla and face (green) ancilla. These can be used to perform 
local stabilizer checks or to mediate many body interactions between edges from 
physical 2-local interactions as described in §4.3. 



does not change the topology of the underlying manifold. Thus, we next add a kinetic 
energy term to the potential, splitting the degeneracy of Hg and reducing to a final 
groundstate capturing homology. 

For each face /, the face Hamiltonian Hf is defined as follows. Orient / according 
to the orientation of the manifold underlying F. Label edges by df = J2k=i ^kSk for 
Ok £ {1, ~ !}• Then we define 



(5) 



With these choices, = for all faces / and vertices v. For the two edges 

incident on a given vertex will be in the boundary of some face, and after correcting for 
orientation conventions this commutativity check reduces to Z^Z'"^] = 0. (See 

Figure f.) Hence, for some constant > 0, we might define iJxE = ^S/g.^^/- 
to commutativity, the kinetic energy Hamiltonian respects the groundstate degeneracy 
of Hq. Label H = Hq + i?KE- We next show that the dimension of the groundstate 
degeneracy of total Hamiltonian 

H = Hq + Hke (6) 
(over C) corresponds to the number of elements of Hi{r,¥d). 

3.2. Homology class groundstates 

The goal of this section is to associate the degeneracy (dimension) of this groundstate 
oi H = Hq + -f^KE to #i?i(r,F(i). We accomplish this in two distinct cases for the 
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manifold underlying T: 

(i) The manifold is orientable, compact, and has no boimdary, so that H2{T,¥ci) = 

(ii) The manifold is compact with boundary and has H2{T,¥d) = 0. For homology is 
a homotopy invariant, and such a surface retracts into its one skeleton. 

Assertion: Let "Hioop denote the groundstate of Hq and H[^] ~ ©,,e[Lj]C |?7). 

^loop = ffiwekcraC|tj) = ©H6Hi(r,Frf) (7) 

Throughout this section, let tt = #G^^ SgGG S- Suppose either Case i or Case ii. 
Then for each [cu], the restriction of tt to H-Ilj] is a rank one projector whose (nonzero) 
image is an element of ker {Hg + Hxe) = ker H. 

To verify this, suppose \lu) is the computational basis state of some cycle w £ Ci 
(i.e. doj = 0.) Then we may also speak of [uj] E Hi{T,¥d), is in the groundstate 
of Hq. Label 

= TT/k) = (#G)-i^.gH. (8) 

It suffices for the Assertion to show the following. 

• li LUi and UJ2 each lie in [lo], then \[uJi]) and |[a;2]) differ by a global phase. 

• If \u) ^ 0, then \[lu]) ^ 0. 

This suffices to see the restriction of tt is a rank one projector, since the first item 
demands the rank < 1 and the second demands the rank > 1. 

We begin with the first item, writing cui — uj2 = E im 82 ■ Since the underlying 
manifold of F is orientable, suppose for convenience that all faces / have positive 
orientation. Then for 77 = Y.fes{v) f P'^* 9v = UfeSM 9f, implying \uji) = \uj2). 
Note that 5,,7r/ tt/^,, = tt/. Thus \ [uji]) = tt/^,, \u}2) = tt/ |w2) = \ [i^2])- 

We next demonstrate that tt/I-Hj^j has rank > 1. As discussed in §2, it suffices 
to show that the trace of this projection, when restricted to the subspace 
which it preserves, is nonzero, and that immediately follows if (f'ld^ ^ demands 
^ = 1. For all other elements of V{n,d) are traceless when restricted to T^ioopi since 
gri\u!) = \u} + dri). Case i and Case ii differ somewhat. In each case, multiples of 
the identity in G/ are products g,, for [77] e H2(T,¥d)- In Case i, besides the empty 
product of the g/ we also produce multiples of /d" as the full product O/ej^ff/' 
< fc < d — 1. This corresponds to H2(T,¥d) = ¥d- Yet for these products = 1, 
as may be verified at an individual edge. In Case ii, there is no nontrivial product 
of the gf which produces a multiple of the identity. This is due to the retraction 
demanding H2{T,¥d) = 0, the second homology of the one complex we may retract 
onto. Colloquially, taking a sum of all faces will force a boundary edge to be acted 
on nontrivially by 5/ for the single face it bounds. Thus in Case ii the only multiple 
of the identity is the trivial product of the g/, and ^ ^ I tautologically. In each case, 
TTflui^] is not traceless and hence has rank at least one. Given the last paragraph, the 
rank is exactly one. 

Retracing the argument above, we may compute the image under tt of the code 
space of Gv is ©[(^]g//j(r,F<i)'C which is also the code space of G. Since tt/ is a 
rank-one projector when restricted to each Ti-iui] ■, we have the following. 

dime ( groundstate of i?) = #Hi{V,¥d) (9) 
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3. 3. Groundstates on a punctured disk 

In practice, constructing physical realizations of Hamiltonians corresponding to two- 
complexes without boundary is daunting. It is possible to simply identify opposite 
qudits on the square fundamental domain of 5*^x5'^= /I? ^ but this would require 
some sort of nonlocal coupling on the boundary in addition to the standard lattice 
coupling. Given a lattice Hamiltonian that arises from electromagnetic coupling, one 
could speculate about some kind of apparatus (perhaps involving fiber-optic cabling 
[2f]) which allows for interactions between boundary qudits. 

Alternately, we might modify the homological groundstates to allow for a surface 
with a boundary curve and punctures. Consider a cellulation F of a disk with k 
punctures. An example with A: = 2 is shown in Fig. 2. Label the j-th puncture face 
/j which has the same orientation as F. Also label the outer boundary of the disk dT 
and the boundaries of the j-ih puncture Analogous to the previous construction, 
the Hamiltonian on F is defined H' = Hq + H'j^^. Here the kinetic term is modified 
so that the set of face operators does not include operators on the punctured faces 
/j, i.e. H'^^ = hJ2f(rjr, Hf, where J^' = !F \ {f £ Ll'-^ifj}. Consequently, there are 
edges on the boundaries dfj that are acted on by X operators from faces on one side 
only. Another way to see this is that all edges of the dual cellulation that cross the 
boundary 9/j share a common vertex located at /j in F. Each edge in F has two 
vertices in V, hence the product over all vertex operators is: 

Y[9v=Id-- (10) 

Not every edge in F borders two faces in T', however, and the product over all face 
operators is: 

k 

l[ gf = Car{X)l[C,{X), (11) 

feJ^' j=i 

where CgriX) = Ylej£dr^°j ^^"^ Q(-^) = Yiejedf -^"j ■ "^^^ orientation Oj = 1 
if the edge ej is oriented in the same direction as the boundary on which the edge 
resides, and = d — 1 if the orientations are opposite. 

First we argue that the code space in nonempty. Recall, the code states are 
defined as -1-1 eigenstates of the stabihzer group G' = {{gf\f £ LI {gv})- The 
operators {{gv}) and ({.g/} commute and the only additional relations obtained from 
the stabilizer group, embedded in Eqs. 10,11, guarantee that {G' f]£,^ Id") = Id"- We 
next show that the code space is Hgr ~TC{k,d) by considering the action of operators 
that commute with any member of G" but act non trivially on Tigr- One such set of 
operators arc non trivial ¥d valued cycles on F generated by {Cj{X)}. We do not 
include the non trivial cycles generated by Ggr{X) because by Eq. 11 their action on 
the code subspace is not independent but can be generated by the cycles around the 
boundaries of the punctures. A non trivial cycle on F is generated by a string of Z 
operations along a path Path{j) that begins on an edge of dfj and ends on an edge 
of dT without touching other edges on puncture boundaries. We denote the generator 
of a such a cycle Cj{Z) ~ Y\ejePath{j) '^°] where Ofc = 1 at the edge € 9/j if 
eu and 9/j share the same orientation and = d — \ otherwise. The other Oj are 
chosen in a consistent way such that C(Z) ji\ =0. The operators on cycles satisfy 
the commutation relations G(Z)°:G{Xfj = £,''^Gj{XfCj{Z)°' for a,6 G as is easily 
verified by considering the action on the one intersecting edge e S 9/'. As such the set 
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Figure 2. An oriented two complex F, which is a cellulation of a two punctured 
disk encoding two logical qudits in n physical qudits. Vertex operators are 
k local where k is the valence of the vertex whereas all face operators Hf are 4 
local in this example. Ground states are +1 eigenstates of the stabilizer group G', 
but not all the stabilizer generators are independent. There are two independent 
non trivial cycles on T which can be generated by closed loops of X operators 
around the boundaries df^ and 9/^. Similarly, there are two independent non 
trivial cycles on the dual F which can be generated by strings of Z operators that 
connect two independent pairs of boundaries of the complex. Shown are the Pauli 
group operations Z^X^ on qudit a and Z'^X'~ on qudit b. 



Rj — {Cj{Z)"-Cj{X)'')}'^ jjiLp generates a representation of the Pauli group V{1, d). For 
sufficiently spaced punctures, all paths Path{j) exist and the group R = {{Rj}'j=i) 
forms a representation oi'P{k,d). We then find that the ground subspace of H' encodes 
k qudits and the set Rj performs local Pauli group operations on the j-th qudit. 

In a lattice implementation of our model Hamiltonian, the punctures may arise 
as physical defects in the system. Coding operations that correspond to cycles around 
defects vividly illustrate the fact that even short ranged correlators (short relative to 
the system size) in a topologically ordered state can have non-trivial values. 

4. Quantum Memory: Input/Output and Error Detection 

We next describe how one might exploit abelian anyons as quantum memories; the 
qubit case has been studied thoroughly [7]. In the new setting of prime level qudits, 
we must treat storage and retrieval of quantum data. It is also possible to generalize 
earlier discussions of stabilizer operations on topologically stored data while in code, 
but we will not treat that topic here. 
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Storing qudits 

Placing quantum data into such a |[a;]) is difficult. For large lattices, this would 
be a special case of the qudit state-synthesis problem. Universal circuits of two-qudit 
operators capable of reaching arbitrary n qudit states are known to scale exponentially 
with the number of qudits [4]. In this section, we propose an alternative which requires 
a number of stabilizer measurements that is linear in the size of the lattice and also a 
sublinear number of entangling gates. 

For an oricntable, connected, compact surface of genus g, it is well known that 
Hi{T,¥d) = (Fd)^^. (See e.g. [15].) We next describe how one might transfer a qudit 
\tp) stored within an ancilliary copy of 7i(l, d) to the topologically ordered groundstate 
eigenspace of H, say Hgr = " . 

The suggestion for encoding is as follows. We begin with \tp) ~ J2j=o '^i b)- 
Choose a copy of C Hi{T,¥d), and let [cu] correspond to 1 G F^. Choose lu G [uj], 
preferably with as few nonzero (excited) edges as possible. Now joj is also a cycle for 
< J < d — 1, and by our choice {[ji-u] = ji^llj^Q contains distinct homology classes. 
Using whatever unitaries are convenient, we form 

d-l 

E 



^) = Va.ljc.) (12) 



For example, on a toric F one might have n-sites and choose a vertical or horizontal 
cycle on a square fundamental domain. Then the appropriate unitary would cost 
0{y/n) gates. Our goal is to construct [^storage) = X]j=o '^i Ib^D- For the remainder 
of the construction, note that all intermediate states are in the code space of the 
stabilizer Gv = {{gv}) Q 'Pin, d). Hence, we may correct for errors in this code at any 
time. Also, the scheme below might be thought of as arising from an error correction 
to the stabilizer G/ = ({.g/}). Nonetheless, only [^'storage) is in the code space of the 
full stabilizer G. Arbitrary local errors are correctiblc in the code space of G since the 
normalizcr of G contains {Ze,Xe ; e e £} [!), 11]. Since this is clearly false for Gy, 
one should perform the initialization above as quickly as possible. 

We suppose an ordering of the faces f G say /i , /2 , . . . , Jl, such that for each 
fixed £ the boundary of fi contains some edge eg which (i) is not within the boundaries 
of /i, f2, ■ • ■ , fi-i and (ii) does not intersect the support of u. This is not possible 
for the last face /l, but we only require this condition for 1 < ^ < L — 1. To store the 



qudit beginning with 



'0 ) , we apply the following steps for each fg 



• Measure the eigenvalue of g/^, e.g. using an ancillary qudit. (See §2.2.) The 
eigenvalue A will be an element of {^^}'^Zq. 

• If A = 1, then the state has collapsed onto the stabilizer {{gv} LI {gf^ ; 1 < 
k < £}) (by induction.) Else, measuring accidentally performed the collapse 
Pj = (l/t^) X]fc=o J which is in fact a projection | Let be the isolated 
edge as above. Since Z^^gf^ = ^''gf^Z^, we see that Z^^Pj ~ PqZ^^. Thus an 
appropriate power of Zf,^ will fix the projection onto the unwanted eigenvalue so 
that the final state lies within the the +1 eigenspace of ff/^. 

Applying the process of the last paragraph clearly produces an element of Hgr ■ 
The applications of Hf, Pj, and also Ze all respect Hyi^] for < j < d — 1. Note 

t Why is this a projection? Consider the unitary h = gf,, and consider projection onto the stabilizer 
of {{h}). 
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that Hgr n Ti-y^] = C If 5* denotes the superoperator of the above sequence 

of measurements and unitary maps, then equivalently we have shown S'(7iioop) ^ 
Ti.gr ri Ti-ijuj]- Equahty is immediate after noting S \ = \ 

However, the effect of the superoperator on relative phases is still unclear. Given 
the global phase on jw), there is a natural global phase on |[a;]) — tt Iw). With the 
argument above, we have actually verified that S \0) ^ c"^" |[0]), S\uj) = e'"^^ 
S \2uj) = e"''^ |[2u;]), etc. Thus perhaps jV'storage) = Sj=o e^^^ctj We argue that 

all of these relative phases are in fact equal. For in terms of the observed eigenvalues, 

L-l L-l L-l 

s = \{zt^'P,{h) - X{P^{fi)zi' = ^X{zt' (13) 

3 = 1 j=l 1=1 

By choice of the support of \lo), also the support of \iuj), we have Jlfci^ ^1^"' li"^) = 1- 



Thus, the applying the superoperator S to i/y produces jV'storage) = X]j=o '^i lli^l)' 
given that we may choose the {ei}^~^ to be disjoint from the support of lu. 



d-1 



4-2. Retrieval 

Thus we next consider retrieval of a qudit stored as in the last subsection, i.e. swapping 
the data in a topological qudit with that encoded in some ancilla qudit. Physically, 
this is more intricate than encoding, which amounts to creating a cycle class \uj) and 
then applying stabilizer corrections for generating Gf. 

For retrieval, the central point is that we may apply a logical X operation to the 
encoded qudit using 0{^/n) gates. To sec this, for uj = '^^rLf.e let X®'^ = (8)egfX"=. 
This might be thought of as a creation operator of an excitation of the loop 
and moreover X®'^ is an element of the centralizer of G not contained within 
G. As such, it preserves the code space, and one readily verifies that it must 
map \[joj\) f-^ [[(j + l)^]), up to global phase. Hence, we may apply controUed- 
X operations targetting the topological qudit using 0{-/n) physical controUed-X 
operations. 

We next consider a controUed-X operation controlled on the topological qudit and 
targetting an ancilla. One must choose a cycle in the dual complex to F according to 
w, say 77. For example, an earlier work [7, Fig. 3] depicts a picket fence dual to a loop 
generator of the first homology group of a torus. In order to perform the required 
controlled-X, follow these steps. 

• Prepare a second ancilla. Then prepare this second ancilla so that the Z eigenstate 
of the ancilla measures Z'^^ . 

• Perform the controUed-X contingent on this second ancilla. 

• Disentangle, i.e. reverse the qudit gates of the first step. 

Consquently, we can perform either controUed-X to or from the topologically encoded 
qudit. 

The ability to perform a two-qudit controUed-X gate implies the ability to 
perform controlled modular addition. The composition begins with a single controlled 
increment triggering when the control carries |1), continues with two controlled 
increments when the control carries |2), etc. The entire circuit thus realizes a controlled 
modular addition in a number of controlled-X gates roughly the triangular number of 
d. Controlled modular subtraction is similar. 
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Finally, modular addition and subtraction allow us to SWAP the topological qudit 
to an ancilla. For bits, the standard three CNOT swap relies on the fact that CNOT 
exclusive-or's one bit to another. Thus the CNOTs perform 6162 i-^ bi{bi © 62) 
^2(^1 © ^2) ^ ^2^1- In like manner, we may perform suitably controlled and targetted 
additions and subtractions for the following sequence of dit operations: 

did2 di{di+d2) ^ i-d2){di+d2) ^ {-d2)di (14) 

Hence, modifying gates so that a control symbol with a + or — target means to add 
or subtract the control respectively, we have the following diagram: 









-1 









We have not described how to complete the gate \j) ^ \d — j) on the topologically 
ordered state. Rather than do so, we claim the top line as the ancilla. This is also 
improves the cost of the controlled additions. 

4-3. Modified constructions using ancillary qudits 

In the quantum circuit model of computation ancillary particles are often used as a 
means to assist in gate operations and as an entropy dump during error correction 
cycles. In the context of surface codes it is tempting to borrow this idea and place 
qudits at the center of each face and on each vertex of the cellulation V, so that the 
appropriate stabilizer checks might be done in place (see Figure lb). Recall, any 
state may be projected into the groundstate of the topologically ordered Hamiltonian 
H = Hg+Hxj?^ using stabilizer checks to the Pauli tensors {gv}U{gf} C Vin, d) (§3.2.) 
Each individual stabilizer check may then be performed using a certain sequence 
of two-qudit gates and a neighboring ancilla (§2.2.) In fact, this basic observation 
presents an auxilliary Hamiltonian which also computes the same topological order 
as the original. Namely, on the face-edge-vertex qudit system, one may build a 
Hamiltonian which is in the groundstate iff all the stabilizer checks g^ and gf are 
satisfied. For suppose we use for the sum gate targetting the qudit of vertex v 
and take rii, to be the qudit number operator on v. Then 

H. = n n (Se)''n„ n (^er' n (is) 

[*,■(;]— e [v.^] — e [*.,v] — e [v.^] — e 

Then \ip) is in the groundstate of Hy iff gy [ip) ~ lip). Similarly, fix a face f € J- with 

9f = EUi "j^J fo"" ^ '^-l}- We take Ff = UUi(^d)e, and Uf = n'=i(S4)"' 
for E{ the sum gate targetting the / qudit. Then for n/ the number operator of the 
qudit at the center of the face /, we label 

As before, we see that \tjj) is in the groundstate of Hf iff 5/ = Thus for /i > 
and U > 0, if = UY^^ Hi, + h^j: Hf, then the groundstate of H is also the code 
space of G = {{gv} U {g/}), i-c the topologically ordered groundstate spanned by 

{|H) ; H ei/i(r,F,)}. 

We finish this section describing another utility for ancillary particles, namely 
to mediate many body interactions present in the Hamiltonian H (Eq. 6) using 
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more physically motivated binary interactions. Consider the vertex constraint term 
Hy ~ ^{dv + gl) where the valence at that vertex is k. This fc- local interaction 
can be obtained as a perturbative limit of 2-local interactions between each d-level 
qudit incident at v and a fc-level ancillary qudit a located at the vertex. Begin with 
a local Hamiltonian for the ancilla Ha = — E'a 10)^^ (0|, and a perturbing interaction 

Va = Jv Ylt=A^er <8) k — 1) (''I + h.c), where ^ \Jv\ and the edge orientations 
give Oj = 1 if Cj = [*,v] and Oj = d — 1 if Cj = [v,*]. By construction, the 
lowest nontrivial, i.e. non identity, contribution to coupling in the ground subspace 
T^gr = \0)a (0| (^0 + V) |0)^ (0| is the effective Hamihonian Hy^s = U{Hy + 0(e)) 
where U = {—l)^Ea{Jv/Ea)'' with an error term of norm ||e|| ^ 1. By judicious 
choice of sign(Jt,) it is possible to fix J7 > 0. A similar argument applies to building 
the face constraint Hf using a j-level ancilla b located at face / to mediate interactions 
between all j edges on the boundary of /. Here we choose Hi, = — £'{, |0){, (0| and 
Vb = ^/Er=i(^e; ® \r-l){r\ + h.c.) such that iJ/cff = h{Hf + 0(e)), where 
h = {—iyEb{Jf/Eby. These mediator qudits could be placed on all the vertexes 
and faces of T to build an effective Hamiltonian in the subspace spanned by states 
with all ancillae in the |0) state. 

An argument in Ref. [ 1 9] suggests that an effective Hamiltonian between spins on 
a two complex can be built using such mediating interactions that closely approximates 
a target Hamiltonian projected to its ground subspace Hgr = PgrHPgr. In the present 
context this would imply that for sufficiently large energies Ea , Et both the degeneracy 
of the ground subspace oi H — Hg + i?KE as well as the energy gap to the excited 
states could be accurately approximated by a model built from a sum of effective 
vertex and face operators. An analysis regarding the validity of such constructions for 
topologically ordered states is wanting, but is outside the scope of this work. 

5. Other Homological Groundstates 

We have originally presented the case of groundstates for Hi{r,¥d) for d prime, in 
order to present the new orientation conventions in the simplest possible context. This 
section describes a construction for homological order on dits whose number of levels 
is not prime but rather a prime power. Homological order for arbitrary composite d 
follows immediately through a tensor product of the prime-power Hamiltonians. 

5.1. Homology Fj^i stabilizer codes 

The hypothesis in the main text has been that qudits have d levels, for d a prime so 
that each \j) is associated to an element of F^. Recent work [11] extends stabilizer 
techniques to the finite fields of order d^, i.e. ¥^e, which exist for any £> 1. 

The generic F^j* constitute all fields F with #F < oo, so this is (perhaps) 
the most general field for which a stabilizer code makes sense. The most typical 
construction of is to consider the polynomial ring ¥d[x] and divide out relations 
in the ideal generated by some irreducible polynomial f{x) ~ + ag^ix^^^ + . . . + aq, 
Qj E ¥d. It is typical to label a E F^e as the adjoined root corresponding to 
the class of x. The Galois group of the extension F^e over F^, say K, then acts 
as permutations of the roots of f{x). Note that ¥^1 is a vector space over the 
scalars F^. Moreover, multiplication by any fixed a G F^f may be viewed as a F^- 
linear map, with an associated matrix with entries in F^. Computing the trace of 
this matrix creates a map Tracep^^ : ¥j^e ¥d. Another characterization is that 
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Tracep /F^^ (x) = X^/teifl'^ ' To ground the discussion, let us review not extensions 
over finite fields but rather Tracec/H(^;) ~ z + 'z ^ 2Re{z). The complex conjugate 
is the Galois action that interchanges i <-*■ —i, for C = M[a::]/(a;^ + !)■ We might 
instead form a 2 x 2 matrix for multiplication hy z — x + iy, which results in 
=x\0) (0| -y|l) (0| -y\0) (1| +x\l) (1| with trace 2x. 
For Wj^e extending the Galois group K is cyclic of order £, generated by si-^x*^ 
for X £ Fjjf . Now K generates a one-qud^it unitary J/^ by U/^lx) = \k ■ x), and the 
corresponding diagonal unitary on the entire lattice will be denoted Uk- 



5.1.1. Fourier transforms for F^jf Having reviewed the machinery of finite fields, we 
next review what one would mean by a stabilizer code of Pauli matrices indexed by it 
[11]. Since our earlier qudit operators X and Z for ¥d had order d, we might instead 
claim to have constructed an X operator and a Z operator for each a G Fd, i.e. X° 
and Z'^ . For Fj^?. we do not take operator powers. Label 7i(l,c?^) = ©aeF^^C{|a)}. 
Then suitable definitions are as follows, where we define ^ = exp(27ri/(i). 

/ X{a)\h) - \a + b) 

For ^ = 1, this generalizes the powers of earlier Pauli operators. Furthermore, with 
these conventions we have a commutator relation: 

X{a)Z{b) = i^''^''''<i^''''^'^''^'^Z{b)X{a) (18) 

Finally, let n{n,d^) = H{l,d'^)'^". In a slight abuse of notation, for a,b E (F^f)" we 
will write a»b = Trace]F^j/F^(ao6o + ^i^i + • • • + a„_i6„_i). Then we may generalize 
the earlier commutatator formula for Pauli tensors as 

[Z{bQ) ® Z{bi) ® ■ ■ ■ ® Z{bn-i)] 

[X{aQ)®X{ai)®---®X{ar,-i)] = , . 

e"'[^(ao) ® X{ai) ®---® X(a„_i)] ^ ^' 

[Zibo) Zibi) ■ ■ ■ ® Z{bn-l)] 

Given this relation, one may define ^^(n, d^) to be that group generated by products of 
Pauli tensors indexed by F^t , as above. Continuing, we may consider sets of particular 
Pauli tensors X(ao)2'(&o)<8i- • ■®X{an-i)Z{bn-i) and consider the stabilizer subspaces 
of the subgroup G C 7^(n, d^) they generate. The error lengths of such code are studied 
in detail [11]. 

Before considering which of these stabilizer codes arise as topological orders, 
we add a point omitted in the original treatments. Namely, we wish to propose 
quantum circuits for the appropriate stabilizer checks. We suppose the existence 
of a number operator measurement which can output classical values in the finite 
field, say abusively n = X^ogf <> ^ I'') ("^l- Then as with F^, stabilizer checks would 
follow given an appropriate Fourier transform T^i : TC{n,d^) ~> T-l{n^d^) which 
maps X{a) eigenstatcs to \a). This leads one to guess we should define Tj^i \a) = 
id')-'/'Etew,,Zia)\b),i.e. 

Ts = {d')-^'^ E '^'^'"'^ \h) {a\ (20) 

However, note that X{a) now has degenerate eigenspaces when i>2. Thus, it is not 
clear the the above equation actually defines a unitary matrix. 
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We briefly comment on why unitarity holds. For convenience, let us drop the 
subscript from the appropriate trace maps. A computation reveals that the unitarity 
assertion is equivalent to knowing that for any fixed a S which is nonzero, 

^Trace(a6) JL q (^21) 

Since a 7^ has a multiplicative inverse, this amounts to 

^Trace(&) ^ g ^22) 



E 



Now suppose we use a to denote the formally adjoined root of f{x) in ¥^ii = 
¥ d[x\ / {f {x)) . Then since every equivalence class may be written as a polynomial 
of degree less than we see that {a-' j^zj is a basis of F^* over F^. In terms 
of the last basis, we might express a generic polynomial class in coordinates as 
h = ht-\X^~^ + hi-2X^~'^ + • • • + 60 for hj € F^. Then Equation 22 becomes 

j^Trace(a*-i)j()f_ij^Trace(Q'-2)jhf_i . . . ^^Trace(l) jbo (23) 

This will in fact be zero, unless all Tracc(a-') ~ mod p, 0<j<^— 1. A standard 
construction in field extensions is to form the discriminant of a basis, for our basis 
A = X^^^Lo Trace(Q;-''''''') |fc) For a given basis, it is not possible that this matrix 
A have determinant zero in F^ [17, Thm2.37,pg.61]. Since the first column of A can 
not then be zero, all Trace(a'') may not be zero, and unitarity of T^it follows. 



5.1.2. Homological order for The chain complex for computing HiiV ,¥^1) 
extends our early discussion by allowing for coefficients of the vertices, edges, and face 
to be within ¥j_i , which in context is (. copies of since only the additive structure 
is relevant. Yet the previous section has nontrivally extended our definition of X and 
Z operators to account for field multiplication, and these operators may be used to 
form a homological order on the physical system in which qud^its (with d^ levels) are 
associated to the edges of F: 

• For each vertex, we may again set g^ = ne=[*,t,] ^e(l) He^l^,*] ^e(-l) and 
Hi, = -{gv + gl)- Then again Ho = UJ2vev 

• Again set gf = (oi)Xe2 (o2)Xe3 (03) . . . ^^^(op), where df = X^Li o^^^ • 

Hf = — (3/ + .9y). Given the generalization of the commutators of the new X and 
Z operators, [Hf, Hy] = for any f,v. Then for h > 0, -ffxE — ^^fej^^f. 

• So if = Hg + i?KE- A similar argument to that given before produces a basis 
\[uj]) of the groundspace of H, as [lo] runs over all elements of fl"i(F,F(j«). 

• These groundstates may again be viewed as a stabilizer code of G = ({ffu, <7/}) ^ 
P{n,d^). Stabilizer checks can be performed as before (see §2.2). The only 
required modifications are that the quantum circuit uses the new Fourier 
transform over F^£ to measure X{a) operators and powers thereof and the number 
operator measurement now takes values in F^f . 

We close with one further comment. Recall [/„ which act on each qud^it as 
J7k |a) = \k ■ a) for k the generator of the cyclic Galois group of ¥^1 extending F^. 
Now UkH — HUk,, as one can verify directly using Hf and Hg. Thus we may view 
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or more generally the Galois action as a symmetry of the topologically ordered 
groundstate. Also, tt^ = ^~^S^=o^k '^i^l ^^len act as a projection collapsing the 
groundstate associated to elements of Hi{r,¥^e) onto the groundstate parametrized 
by Hi{T,¥d) as constructed in §3. In terms of Hamiltonians, tt^ projects onto the 
groundstate of H,^ = —{11^ + Ul), whose physical significance is unclear. 

6. Z/dZ Gauge Theory and Anyonic Excitations 

In our treatment of code subspaces, we have used the isomorphism between spins on a 
surface and one-chains on a two complex to label the ground states of the Hamiltonian 
H in terms of homology equivalence classes. The language of cell complexes also carries 
over to describe the excited states. If we identify the ground subspace of H as the 
vacuum then excited states are labeled by valued boundaries of one chains on the 
complex r or the dual complex F. These excitations can be viewed as massive particles 
with definite statistics. 

In this section we show by construction that our model is a Z/dZ gauge theory 
with quasi-particles corresponding to dyonic combinations of charge and flux. These 
quasi-particles have abelian anyonic statistics. We provide an algorithm in terms of 
an interferometer circuit for measuring components of the scattering matrix. 

6.1. Stabilizer errors as abelian anyons 

Consider a two-complex T with a physical system of qudits associated to each edge 
and a topologically ordered Hamiltonian H as above. We have already seen how 
to associate a basis of the groundstate eigenspace with elements of Hi{r,¥j). As 
stabilizer states, it is well known that the groundstates are entangled. Abelian anyons 
arise as entangled excitations of this system. In the qubit case, such excitations always 
arise in pairs [12]. In our generalization, this is also true, and an excitation \ j) is always 
paired to an excitation \d ^ j). 

The linear algebra for constructing a charge anyon is as follows. First, choose 
two vertices vi and V2 of F on which the anyon should reside with charges j and 
d — j respectively. Choose a chain u; with dcu — jvi + {d — j)v2. Recall from §3.2 the 
projection tt = Tliferdf which projects onto the stabilizer code of all the face 

operators g/ = ®e(^df[Xe)'^ ■ We set 

IV'charge anyon) = 7^ 1^) (24) 

The resulting state is an excited state of Hq whose eigenenergy is 4[/(l - Re(e2*'^^'/'*)) 
above ground. It is not independent of the choice of w, and this in fact allows for 
an interesting geometric interpretation of the error length of the associated stabilizer 
code [7]. 

For let oji and 0J2 be two such choices, with 1-01) and \i]j2) the resulting anyon 
states. Then loi — 0^2 is a cycle, and 

!</>) ^=''^(K)-M) = |^i>-|V'2) (25) 

is the ground state eigenket associated to [011 — ^2] G i?i(r, F^). Hence, if we encounter 
such a charge anyon excitation which has sullied a qudit encoded in the groundstate 
of H , then correcting it amounts to choosing an cancelling anyon or equivalently to 
choosing a cycle on F. If the dual charges of the anyon are separated by roughly 
half the diameter of the two-complex, then this choice is likely to cause an error. Yet 
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for nearby dual charges one might reasonably guess [lui — 1^2] — [0]. In particular, if 
r were to cellulate the square fundamental domain of a torus using n qudits on the 
edges (implying Q{y/n) qudits on a side,) then we would expect an error length for 
the associated stabilizer code to be roughly 0{y/ri) [7]. 

Similar comments apply not only to charge anyons but also flux anyons [12]. Here, 
one chooses a path in the dual complex to F, i.e. a sequence of connected faces. Let 
|[0]) be the homological groundstate associated to [0] € ifi(r,Fd). A flux charge of 
multiplicity j on the endpoints of the face path is associated to 

I V'flux anyon) = ^vdf;' 9%^ ' ' ' 5// I M) (26) 

where 7r„ = i#V)~^ J2v£V 9'" ^^'^ path consists of faces /i,/2,...,/£ with the 
signs allowing for orientation. The flux anyon theory follows quickly by considering 
the charge anyons of the dual two-complex to F, say F. Faces of F become vertices 
of F while vertices become faces, and the graph of F arises by connecting vertices 
corresponding to incident faces of F. Suitable hypotheses on the cellulation of the 
underlying two-manifold of F will cause this dualization procedure to be well behaved 
[15], and one might associate charge-anyonic observation of flux anyons and vice versa 
with pairings exploited in the proof of Poincare duality. 

6.2. Quasi-particle statistics 

We next wish to study such anyon states, i.e. errors of the stabilizer code as 
above. New notation for the excitations follows. A charge a G Z/dZ at vertex v 
is labeled by the state |(a, 0; (w, — ))) such that ((a, 0; (w, — ))| g„ |(a, 0; (w, — ))) = 
Similarly, flux b G Z/dZ at face / is labeled by the state |(0, 6; (—,/))) such that 
((0,6; (— ,/))|5j 1(0,6; (—,/))) = A dyon refers to a bound state of charge and 
flux at vertex v and face / neighboring each other, i.e. [v,*] G df or [*,v] G df 
and (a, 6) G (Z/dZ)^. The state of such a dyon in Hilbert space will be denoted 
I (a, 6; (v,f))). For simplicity we restrict our discussion to simply connected compact 
surfaces with boundary such that the ground (vacuum) state is nondegenerate § 

Pauli-group elements local to a single edge of F produce dyons of the topological 
order in particle anti-particle pairs. To see this, note that the operator acting at 
edge e = [wi, V2] creates a pair of boundaries on the vertices, one with charge a at vi 
and other with charge d — a at V2. We name the charge d — a particle an anti-charge 
to a. Similarly, the operator creates quasi-particles located on the two faces /i and 
/2 that share the edge e on their boundaries. Let face /i be the face with opposite 
orientation to e. Then the flux at /i is b and the anti-flux at /2 has the value d — b. 
A product operator X°'Z~^ acting on edge e creates the dyon (a, b) with charge a at 
vertex vi and flux b at face /i (sec Fig. 3a). When it might be clear from context, 
we will drop the particle location labels (f,/), e.g. particle anti-particle pairs might 
be written as |(a, 6); (—a, —6)). The mass of a dyon is given by the expectation value: 
ma,b = {ia,b)\H\{a,b))- Eq = 2C/(1 - Re[C"]) + 2h{l - Re[e^]), where Eq is the 
vacuum energy. The energy to create a particle antiparticle pair is twice this value. 

Prior work in continuum field theory has considered dyon excitations in which 
charges and fluxes take values in Z/dZ. The interactions described by a Z/dZ gauge 

§ In general |(a, b); {v, /)) describes an equivalence class of pure states which results from applying 
X^Z^^ to any groundstate. For a degenerate vacuum, particle creation, followed by braiding and 
annihilation can result in non trivial logical operations on the code subspace. 
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theory are completely characterized by the foUowmg rules [211]. 

|(a, b; {v, /))) X |(a', b'; {v, /))) ^\{a + a',b + b' ■ (v, /))) (27) 

7^ |(a, b- {v, /))) \{a, b- {v', /'))) = C' b; {v, /))) |(a, 6; {v' , /'))) (28) 

7^2 |(a, 6; («,/)))|(a', 5'; («',/')))= 

e-''>+''-ma,b;{v,f)))\{a\b';{v'J'))} ^^^^ 

C |(a, 6; («, /))) = !(-«, -6; («,/))) (30) 

r|(a,6;(«,/))) =r''l(a,&;(t',/)))- (31) 

We next review these rules and argue that the dyonic excitations of our Hamiltonian 
satisfy them. 

The first relation is the fusion rule for particles occupying the same location where 
addition is performed modulo d. In the context of our model this rule follows from the 
additivity of boundaries of one chains. Indeed, it is the ability to annihilate particle 
anti-particle pairs by choosing a trivial cycle on F or F that makes correction of local 
errors possible (see Fig. 3b). The next two rules describe the action of the monodromy 
operator TZ which performs a counterclockwise exchange of one particle with another. 
The quantum state of n indistinguishable particles residing on a surface belongs to 
a Hilbert space that transforms as a unitary representation of the braid group -B„. 
If we order the positions of the particles {(wj , /j)}"=i, then the n — 1 generators 
of Bn correspond to the monodromy operator TZ acting on the particle pairs in the 
locations {{vj, fj), (wj+i, fj+i)}"Zi ■ a Z/dZ gauge theory, the irreducible unitary 
representation of Bn is one dimensional, meaning the particles are abelian anyons. 
Notice that the definition of the monodromy operator involves orientation of the path 
taken during particle exchange. For a non oricntablc surface, Z/dZ statistics for > 2 
are not allowed because the clockwise trajectory of particle around another is not 
uniquely defined whereas the phases arc distinguishable except for d~2. 

The braiding of one dyon around another is shown in Fig. 4. Here 
we begin with a state of two dyonic particle anti-particle pairs: j^*) = 
|(a, 6); (— a, — 6)) |(a', 6'); (— a', — 6')) in distinct locations on the surface. The mutual 
statistics are determined by winding one dyon, (a, 6) around the other (a', 6') in 
a counterclockwise sense. This action is described by the square of monodromy 
operator TZ, which exchanges two particles in a counterclockwise sense. A non trivial 
phase is accumulated under the action of TZ^ because the closed loop string operators 
that wind (a, 6) collide with the strings connected the dyon {a',b') with its anti- 
particle. In the example shown in Fig. 4a, the strings intersect at two locations 
where we have the operators Z~''X~''' = ^'^'^X^'^'z'' and X-'^Z'"' = ^^'"Z^^''^-". 
Rewriting these operators with the action of the closed strings (Pauli operators 
with unprimcd powers) first has the advantage that the closed strings act trivially 
provided that there are no other quasi-particles inside the closed loops. Hence we 
have that TZ^ l^*) = ^('^' ''+''' The preceding example illustrated the Aharanov 
phase accumulated when winding one charge around a flux along a trajectory that 
was local, i.e. did not explore the global properties of the surface. Were the flux 
absent, then the path would be homotopic to a point. One can also define this 
phase for trajectories that explore the global properties of the surface, but can be 
continously deformed to a process where one anyon wraps around another. On a 
torus, for example, the following process traces out non trivial cycles for the charges 
and fluxes. Represent the torus as a square with opposite sides identified, labelling the 
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axes of the square Xi,X2. Pick a non trivial cycle along the Xi direction of F, and call 
it Pi . Similarly, pick a non-trivial cycle along the X2 direction of the dual T and call it 
P2. To obtain the exchange statistics, first create the dyonic particle antiparticle pair 
|\[') = |(a, b; {vi, fi)); {—a, —b; {v2, /2))) out of the vacuum state Wind the charge 

a around Pi so that it annihilates with its anticharge partner at site V2- Next wind 
the flux b around P2 so that it annihilates with its antiflux partner at face /2. Create 
another dyonic particle antiparticle pair j^E'') = \ia,b; (v2,/2)); i—o-,—b] with 
particle antiparticle positions reversed relative to Wind charge a around Pi in 
the opposite direction to the first winding so that it annihilates with the anticharge at 
site vi and likewise, flux b along P2 in the opposite direction so that it annihilates with 
the antiflux at face f2- These four trajectories cross at one edge e and the action on 
the state (for one choice of edge orientation) is l^-g) Z'^X'^ZlXl \^g) = \'^g). 
If we embed the torus in R^, then the worldlines described by intersecting strings in 
the above process are equivalent under ambient isotopy to linked world lines on the 
plane which describe winding the charge a around the flux b. 

Identical quasi-particle statistics are determined by exchanging one dyon (a, b) 
counterclockwise with another. Such a process is depicted in Fig. 4b. The action on 
the state = \{a,b)) |(a,6)) can be computed by annihilating particle-antiparticle 
pairs after exchange, creating them again, and comparing the resultant state with the 
initial state \^). We can annihilate the charges on the left side first. Reversing the 
order of the operator that created the dyon there, we have X"'Z~'' = if-^Z^X"' and the 
charges are annihilated by applying . Similarly, the charges on the right side are 
annihilated by applying a string of Z~^ operators. Finally, the fluxes are annihilated 
by applying X" or along the remaining two connected strings. The action on the 
wavefunction is then 7^ | (a, b)) \ (a, b)) = \ (a, b)) | (a, 6)) . 

The particle conjugation operator C in Eq. 30 reverses the sign of all the particles. 
This is realized in our microscropic spin model by reversing the orientation of the all 
the edges on the cellulation. Finally, the operation T in rule 31 rotates the charge 
component of a dyon around its own flux, generating an Aharanov-Bohm phase in 
the process. This is illustrated in Fig. 4b. Here the charge component of the dyon 
(r, s) is wrapped around its flux component in a counterclockwise sense. During this 
operation, there is a collision at the edge where the dyon was created. Rewriting 
the operation on the edge as X'^^Z~'' = ^^^-^r ^-sx'^ gQ ^j^g^^ Xoo^^ operation about 
boundary of the face / acts trivially first, we have that T |r, s(w, /)) = |r, s{v, /)). 

6. 3. Measuring statistical phases 

In any physical construction of a Hamiltonian that admits topologically ordered states 
it will be important to verify the predicted properties. One, albeit crude, observable 
is to measure the energy gap from a ground state to a first excited state. This 
could be done by probing linear response of the ground states to a perturbing field 
that generates local unitary operation at a frequency ujf- For a system with the 
internal Hamiltonian Eq. 6, the expected resonant absorption occurs at frequencies 
lof = 2ma_;,//i. However, as a witness to topological order, this measure is not 
sufficient because there could be another spin Hamiltonian with equal gap that does 
not possess topologically invariant correlations functions. Another more convincing 
probe would be to directly compute the statistical phases in Eq. 29. Operationally, 
this should be done by measuring both the phase (f>j- accumulated when one particle 
(a, b) wraps around another (r, s) and the phase 01 when the particle (a, b) traces out 
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Figure 3. Quasi-particlc excitations on a honeycomb ccllulation. (a) 
Excitations appear in particle anti-particle pairs. Charges(anti-charges) appear 
as boundaries on vertices represented by open(filled) diamonds, and fluxes(anti- 
fluxes) as boundaries on the faces represented by open(filled) squares. The 
total charge and flux of any pair is zero. Shown is a flux pair |(0, c); (0, — c)}, 
charge pairs 0); (i, 0)) , |(fc,0); (— fc,0)) and a bound state of charge and flux 
pairs I (a, b); (— a, — fe)). Notice that strings of the same or different types are 
allowed to intersect, (b) Fusion of quasi-particles. The upper two diagrams 
illustrate corrective procedures to annihilate charge and flux excitations. The 
lower two diagrams illustrate the fusion rules |(i,0)) X |(A;,0)) = |(j + fc,0)) and 
|(0,i)> X |(0,-fc)> = |(0,i-fc)>. 
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Figure 4. Braid relations, (a) Counterclockwise braiding of the dyon (a, 6) 
around the dyon (a', fe'): Tl'^ \{a, b)) \{a' , b' )) |(a, fe)) |(a', fe')>- (b) 

Counterclockwise exchange of identical dyons: TZ \{a, b); (a, b)) = |(a, b); (a, b)). 
In the upper left hand side of the surface is shown the counterclockwise winding 
of the charge component of a dyon (r, s) about its flux component generating an 
Aharanov-Bohm phase according to T |(r, s)) = g*"* |(r, s)). 
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Figure 5. Protocol for measuring quasi-particlc statistics. The green circle 
represents an ancillary particle which performs conditional gate operations on the 
qudit residing on edge e = [1)2, ''o]- The red lines indicate operations which are 
done adiabatically with respect to the energy gap AE. The inset is a simplified 
space-time diagram of the braid. 



the same path in configuration space but does not enclose the particle (r, s). The phase 
difference 0r ~ 0i = '/'top subtracts out dynamical phases and Berry's phases, leaving 
only topological information. We sketch an algorithm for computing this phase using 
operations in accordance with the two complex illustrated in Fig. 5. Adaptation to 
other cellulations is straightforward. 

(i) Beginning from a ground state |5'(0)), prepare a state with two particle anti- 
particle pairs in disjoint regions of the surface: 

|vl'(l)) = \{a,b;{v,J,));{~a,~b;{viJ^))) 
\{r, s; {vo, /o)); {~r, -s; {vi, A))) . 

(ii) Prepare an ancillary qubit a in the state \+x)a = '^^^^'^a + \^)a) ^^^"^ ^^^^ 
qubit to perform the controlled unitary operation /\i{X^^Z^) = |0)^ (0| Id + 

{1\^X~^Z^ (with (r, s) ^ (0, 0)) on the qudit residing on the edge e = [^2, t'o]- 
Measure the ancilla in the x basis and record the result m = ±1. The resultant 
state is \^{2)) = ^(|^'(1)> + |*(1)>), where 

X--ZI\'f{l)) - |(a,5;(«3,/3));(-a,-&;K,/4))) 
s; {v2, /2)); {-r, -s; (wi, /i))) 

is orthogonal to 

(iii) Use a sequence of local spin operations to drag the dyon at location {v2, fi) to the 
location (w5,/5). These operations should be done adiabatically, i.e. they should 
be done using localized control fields with frequency components much smaller 
than the minimum gap energy A_E. In this way no new particles will be created, 
only the component of the wavefunction with the dyon located at (w2j f2) will be 



Qudit surface codes and gauge theory with finite cyclic groups 



24 



changed. Instead of using control fields to perform local spin operations, another 
possibility is to slowly decrease the values of U and h on the vertices and faces in 
the path from {v2, f^) to {v^, /s) so that it is energetically favorable for the dyon 
to follow this path. The resultant state is: |^'(3)) = + (-1)" l^-')), 

where 

l^-') = |(a,6;(t;3,/3));(-a,-&;(t'4,/4))> 
\{'r,S] {V5,f5)); (-r,-s; • 

(iv) Braid the dyon (a, 6; (ua, /s)) in a counterclockwise sense around the location 
(wsi/s) such that it returns to location (wsi/s). The state is now: |^'(4)) ~ 

(v) Perform the inverse of the operations in step iii, again insuring that no new 
quasi-particles are created during the process. The resulting state is; |'I'(5)) = 
-^{1^(1)) + (-l)"^("''+'''')e*xx-''Z| |*(1))), where we have included x, the sum 
of dynamical and Berry's phases that may have accumulated during steps ii-iv. 

(vi) Reprepare the ancilla in the state and perform the controlled unitary 
operation Ai((— Measure the qubit in the x basis. The expectation 
value is: 

{'^a)r = 5(cos(x + 0top) + 

S2rM'>2sM cos(x + (/"top - 2nrs/dfj , 

where 0top — 27r(.sa + rb)/d. 

(vii) Repeat steps i-vi but measure the ancilla in the y basis. The expectation value 
is: 

(0-a>r = i(sin(x + 0top)- 

'^2r,o<^2s,o sin(x + 4>top " 27rrs/d)^ , 

(viii) Perform a similar experiment but this time using a trivial braiding operation, i.e. 
perform the steps in the order (i,ii,iv,iii,v,vi,vii) so that the braid is contractible. 
Then the expectation values are 

Wa)i = ^(cosx + '^2r,o'52s,ocos(x- 27rrs/d)^ , 

(^a>i = ^(sinx- <^2r,o'52s,osin(x- 27rrs/d)y 

(ix) Compute the topological phase (/)top from an ensemble average obtained by 
repeated measurements on identically prepared systems. 

As a simple example, consider the computation of the mutual statistics of charge 
and a flux for d = 2. Setting (r, s) = (0,1) and (a, 6) = (1,0), the expected 
measurement results are {<t'^)t = ^ cos(x + 0top), Wa)T = 0, (cr^)i = ^cosx, 
{<^a)i = 0- If desired, the phase x could be engineered to vary in a controlled manner 
over different trials in order to improve the visibility of the phase shift cf>top- For d > 2 
it is always possible to choose the probe dyon such that d2r,oS2s.o = 0. In this case, (/)top 
is estimated by finding the closest solution to e"^'°p = ((o'a)T+j(o'a)r)/((o'a)i+*(o'a)i)- 
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7. Conclusions 

We have proven the existence of a microscopic spin model that provides for 
topologically protected qudit encodings. This model describes a Z/dZ gauge theory 
with abelian charge/flux dyons as excitations. The construction is quite general, 
allowing for arbitrary cellulations of an orientable surface and encoding qudits with 
any finite number of levels. Suggested adaptations to the standard spin models 
using ancilla for in place stabilizer checks could prove advantageous in any physical 
implementation of such codes. Moreover, with some limited degree of local control 
it is possible to measure the anyonic statistical phases. Given that these properties 
arc difficult to measure in quantum Hall systems, this could provide a novel probe of 
topological order. 
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